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Abstract 


Fluid flows are typically studied by solving the Navier-Stokes equation. One of 
the fundamental assumptions of this equation is Stokes’ hypothesis. This hypoth- 
esis assumes bulk viscosity, 4p, to be identically zero. The Stokes’ hypothesis is a 
reasonable approximation for commonly observed fluid flows; therefore, Navier-Stokes 
equation gives satisfactory results in these situations. However, there are circumstances 
where this hypothesis does not hold good, and therefore, the classical Navier-Stokes 
equation becomes inapt. These situations include the absorption of sound waves, hy- 
personic flows, turbulent flows, and flow of Martian air etc. Reliable analytical and 
computational studies of these flows requires account of bulk viscosity in the governing 
equations. In this article, we provide a brief review of the subject of bulk viscosity. We 
start with a brief background of this topic. Then we discuss the underlying microscopic 
mechanisms that give rise to bulk viscosity effects. It was followed by a review of meth- 
ods available in the literature for estimation of this parameter. Finally, a review of the 


studies that analyze the effects of bulk viscosity in various fluid flows is provided. 
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1 Introduction 


Ever since Sir George Gabriel Stokes (1819-1903) proposed the complete set of equations 
for the dynamics of viscous fluids in 1845 i, bulk viscosity (sometimes also referred as 
volume or dilatational viscosity) has remained as one of the controversial subjects of fluid 
dynamics [2]. The recent surge of interest in Mars missions has made the study of bulk 
viscosity more relevant as the Martian atmosphere consists approximately 96% of carbon 
dioxide gas, which has a reported bulk to shear viscosity ratio of ~ 2000 [4]. The account 
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of bulk viscosity in the analysis has also enabled more accurate modeling of several fluid 
mechanical phenomena BHEN]. On the other hand, in contrast to shear viscosity, which is a 
very well studied transport property, bulk viscosity is still not completely explored. There 
are considerable ambiguities/uncertainties about the nature, effects, and applicability of the 
concept of bulk viscosity. Even for most common fluids, existing experimental values of bulk 
viscosity are spread over a broad range, and widely accepted values are still not available[18}. 
In this article, we provide a brief review of the subject of bulk viscosity. We start with a 
brief background of this subject. Then we discuss the underlying microscopic mechanisms 
that give rise to bulk viscosity effects. It was followed by a review of methods available in 
the literature for estimation of this parameter. Finally, a review of the studies that analyze 
the effects of bulk viscosity in various fluid flows is provided. 

The stress-strain rate relationship (i.e., constitutive relation) for a Newtonian fluid is 


given as follows: 


Ou; ðu 
Oik = =F henna Oik f H (+ E) 
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where, oj, is Cauchy’s stress tensor, 0;, is the Kronecker delta, u; is velocity of the fluid, 


PAn (1) 


x; is spatial coordinate, and the scalar quantity Pihermo is the thermodynamic pressure or 
hydrostatic pressure. The above relation contains two independent coefficients - first is the 
coefficient of shear viscosity (u), which is sometimes also termed as the first coefficient of 
viscosity, and the second is the coefficient of longitudinal viscosity (A), which is also referred 
to as the second coefficient of viscosity. 

The above relation, Eq. (1), can be rearranged as follows by separating isotropic and 


deviatoric parts of the strain-rate tensor: 
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where, u = Gr + A). The coefficient 4, is known as the coefficient of bulk viscosity, and 


it represents the irreversible resistance, over and above the reversible resistance, caused by 
isentropic bulk modulus to change of volume [19}. Its values are expressed in the same units 
as shear viscosity, i.e., Pas or poise. The second law of thermodynamics constrains both the 
shear viscosity and bulk viscosity to have only non-negative values [20]. The values of bulk 
viscosity for common dilute gases at 300 K are listed in Table [1] (21). It should be noted that 
the bulk viscosity of monatomic gases in the dilute gas limit is zero. 

The mechanical pressure, Pech, is defined as the negative average of the diagonal terms 


of stress tensor, as given below: 
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Gas H»/H Gas H»/H 
Carbon monoxide | 0.548 | Dimethylpropane | 3.265 
Nitrogen 0.769 Water vapor 7.36 
n-Pentane 0.896 Hydrogen 28.95 
iso-Pentane 1.057 Chlorine 751.88 
n-Butane 1.13 Fluorine 2329 
iso-Butane 2.00 | Carbon dioxide 3828 


Table 1: Ratio of bulk viscosity to shear viscosity for common gases at 300 K 


Historically, Stokes |1| assumed the bulk viscosity (4) to be identically zero for all fluids. It 
implies that mechanical pressure is always equal to thermodynamic pressure, irrespective of 
the process through which the system is undergoing, i.e., the viscous forces do not depend 
upon the rate of expansion or compression at all. This assumption is known as the Stokes’ 
hypothesis. Later, it became customary to use this hypothesis in fluid mechanics. However, 
Stokes himself did not take this hypothesis as always true. He mentioned that in com- 
monly encountered flows, if analysis with and without considering bulk viscosity produces 
the same results, it would be because of small V-wu rather than u, being zero. Indeed, there 
are certain instances, where bulk viscosity effects are not negligible. The Section [4] briefly 
discusses these scenarios, and some of the past works to model them while considering finite 


bulk viscosity. 


2 Microscopic picture of mechanisms of generation of 
normal stresses in an expanding/contracting dilute 


gas 


To understand the microscopic origin of the normal stress, let’s re-write Eq. [2]for the normal 


stresses acting in the x-direction on yz-plane of the fluid element. 
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The second and third term in the above equation are contributions of shear and bulk viscosity 
to normal stress, respectively. Mechanisms responsible for stresses due to these terms are 


discussed as follows. 


u } utdu‘u+2du 
i—>: — 


Figure 1: Three layers of different velocity in a divergent flow field. Dots represent gas 
molecules and two vertical dashed lines are imaginary boundaries separating these three 
fluid layers. 


2.1 Shear viscosity 


In dilute gases, the mechanism of generation of normal stress due to shear viscosity is anal- 
ogous to the mechanism that produces shear stress, i.e., the transport of momentum in 
adjacent layers of fluid due to thermal diffusion of molecules. Consider a flow field dilating 
only in one direction with three different layers of fluid having velocities, as shown in Fig. 
The left, middle and right layers have a bulk velocity equal to u, u+ du, and u+2du, respec- 
tively, in the direction as shown in Fig.|1| When any molecule from the leftmost layer jumps 
to middle layer due to its random-thermal-motion, it will decrease the average momentum 
of middle layer. As a result, middle layer will be pulled towards the left side. Similarly, if 
any molecule from the rightmost layer jumps to the middle layer, it will increase the average 
momentum of middle layer. As a consequence, there acts a pull on the middle layer towards 
the right side. The combined effect of these two pulls will be a normal stress on the fluid 
element of middle layer along the direction of velocity gradient. Therefore, in contrast to 
what the name ‘shear viscosity’ can mislead to, shear viscosity can not only produce shear 


stresses but also normal stresses. 


2.2 Bulk viscosity 


The pressure in any fluid is summation of two components— 
P= Piinane + Pagal (5) 


The first contribution, i.e., Pkinetic is represents force caused by bombardment of molecules 
on a unit area. The second contribution, Pyiria is caused by inter-molecular forces. It 
depends upon the strength of inter-molecular interaction. It should be noted that the P 
can refer to both mechanical and thermodynamic pressure depending on the reference. The 
mechanical pressure, Pech, is instanteous pressure of the fluid at all situations, i.e., both 
at equilibrium or non-equilibrium. Hydrostatic/thermodynamic pressure, Pinermo, is the 
mechanical pressure of the system under hydrostatic/equilibrium conditions, i.e., when fluid 
is at rest. Therefore, for a system in a state of non-equilibrium, the term ‘thermodynamic 
pressure’ loses its meaning [20]. However, for such a system, it can still be defined as the 
mechanical pressure of the system if the system is brought to equilibrium state adiabatically 
42). 

The origin of bulk viscosity in fluids is due to the fact that the pressure at non-equilibrium 
(i.e., mechanical pressure) is not same as the pressure at equilibrium (i.e., thermodynamic 
pressure). There are primarily two mechanisms responsible for deviation of equilibrium 
pressure from non-equilibrium pressure- (a) The finite rate of relaxation of internal degrees 
of freedom with random translational energy (b) The finite rate of structural rearrangement 
of molecules after a change in the thermodynamic state of the fluid. The bulk viscosity 
due to first mechanism is called as ‘apparent bulk viscosity’, while bulk viscosity due to the 
second mechanism is called as ‘intrinsic bulk viscosity’ [3]. The total bulk viscosity of fluid 


is summation due to the two said contributions. 


2.2.1 Apparent bulk viscosity 


The apparent bulk viscosity appears when the Pkinetic at non-equilibrium is not same as 
Pxinetic at equilibrium due to presence of internal degree of freedoms, i.e., (viz., rotational 
and vibrational). Herzfeld and Rice first suggested that microscopic cause of bulk 
viscosity is the finite rate of exchange of energy between the translational mode and internal 
degrees of freedom. The mechanism can be explained by considering a simple example in 
which a polyatomic dilute gas, say nitrogen, is expanding adiabatically in a piston-cylinder 
arrangement, as shown in Fig. |2| Since the gas is considered as dilute, the virial pressure is 
neglected. 

Let us assume that initially the piston was at rest, and a polyatomic dilute gas was in 


equilibrium at a temperature of 300 K. With the assumption of inactive vibrational modes 


Figure 2: Expansion of gas in a piston-cylinder arrangement 


at the prevailing temperature conditions, the gas has three translational and two rotational 
degrees of freedom. Since the gas is in equilibrium, these all five degrees of freedom possess 
equal amount of energy because of the equipartition law of energy. Now, when the gas 
expands, it does work against the piston and loses its energy. The energy that the gas loses 
comes directly from its translational mode, whereas the energy associated with rotational 
mode remains momentarily unaffected. It causes an imbalance in the equipartition of energy 
among translational and rotational degrees of freedom. At this stage, the system is in a 
state of non-equilibrium, and its current translational kinetic energy is less than that which 
would be if the system is again brought to equilibrium adiabatically. Similarly, its rotational 
kinetic energy is more than that which would be if the system is brought to equilibrium. 
The system tries to regain its state such that the whole kinetic energy is equally distributed 
among all degrees of freedom. In this attempt, the system transfers some of its kinetic energy 
from internal modes to translational modes by means of inter-molecular collisions amongst 
gas molecules. 

The mechanical pressure, Pech, at any point in the fluid is the negative average of normal 
stresses acting on that point. For a dilute gas with negligible virial pressure, it represents 
actual force caused by bombardment of molecules on a unit area, and therefore, depends 
only upon the random translational kinetic energy of the gas molecules. This is also the 
reason why PV work comes at the cost of translational kinetic energy. Kinetic theory of 
gases relates mechanical pressure to translational kinetic energy (Erans) by the following 
relation: oe 

trans 
Pk = 3p (6) 
where VY is volume of the system. Moreover, due to the equipartition law of energy, the 
translational kinetic energy (Frans) of the system at equilibrium is equal to 3/f of the 


total energy (Eyota), where f is number of total degrees of freedom of gas. Thus, the 


thermodynamic pressure can be given as follows: 
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From the above discussion, it can be deduced that for the time duration, when transla- 
tional kinetic energy is less than its value at equilibrium, the mechanical pressure is also less 
than the mechanical pressure at equilibrium (i.e., the thermodynamic pressure). In case of 
compression, the situation would be reverse. Since, the energy supplied to the gas first goes 
to translational mode, instantaneous translational kinetic energy is more than its value at 
equilibrium. Therefore, the mechanical pressure is also higher than the mechanical pressure 


at equilibrium (i.e., the thermodynamic pressure). 


2.2.2 Intrinsic bulk viscosity 


This mechanism was first proposed by Hall in 1948. Consider a compression of dense 
fluid such as liquid water. During compression, two different processes take place simultane- 
ously. The first one is that molecules are brought uniformly closer together (just like zoom 
in/out of a computer graphics). It can be called molecular compression, and it is an almost 
instantaneous process. However, even if the molecular arrangement before the compression 
was a stable one, i.e., the one with minimum intermolecular potential energy, the molecu- 
lar compression alone does not guarantee that the resulting molecular arrangement would 
also be a stable one. The reason for this is the nature of the intermolecular interactions, 
e.g., Lennard—Jones potential or hydrogen bonds. Therefore, a second process also happens 
simultaneously. In this process, the molecules are rearranged or repacked more closely to 
achieve a more stable configuration. Hall identified this process as configurational or struc- 
tural compression. This process involves the breaking of intermolecular bonds (e.g., hydrogen 
bonds) or flow past energy barriers, which stabilizes the equilibrium configuration. This 
is a finite rate process. Thus it is of relaxational nature and is a source of nonequilibrium. 
This mechanism of bulk viscosity is present in all fluids including monatomic gases. Hence, 
monatomic gases at atmospheric conditions have a small (O(10~'°) Pa s) but non-zero 
bulk viscosity. It should also be noted that at hypothetical dilute gas conditions, the bulk 
viscosity of monatomic gases is considered to be absolute zero. 

Monatomic gases in their dilute gas limit do not possess any appreciable inter-molecular 
potential energy, therefore have negligible intrinsic bulk viscosity. Moreover, because of 
monatomic nature, they do not possess any internal degrees of freedom. Therefore, the me- 


chanical pressure at any instant will always be equal to the thermodynamic pressure. Hence, 
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Figure 3: A survey of available methods for estimation of bulk viscosity. 


Tisza's formulation 


Li's formulation 


Kustova's formulation 


we can expect them to exhibit negligible bulk viscosity. Both theories and experiments 
also confirm the same. However, it is possible that factors other than rotation/vibration 
and potential energy, like electronic excitation or chemical reaction can also cause non- 
equilibrium in dilute monatomic gases. For instance, Istomin et al. has shown that the 


bulk viscosity is not zero in electronically excited monatomic gases at temperatures higher 
than 2000 K. 


3 Methods for determination of bulk viscosity 


Unlike shear viscosity, determination of bulk viscosity has always remained a challenging 
task. We present a brief summary of various approaches available for estimation of bulk 
viscosity, 4p, of non-relativistic classical fluids, and QGP and hadronic matter in Secs. 
and respectively, with a particular focus on the former. 


3.1 Determination of bulk viscosity of classical fluids 


A schematic overview of the methods available in the literature for the determination of bulk 


viscosity is shown in Fig. |3} A brief discussion on these methods is given below. 


3.1.1 Theoretical methods 


Theoretically, bulk viscosity of dilute gases can be related to the relaxation time of equili- 


bration processes by Tisza’s |3| formulation, given as follows: 
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where peg is density of gas at equilibrium, a is speed of sound in absence of viscosity, y is 
the ratio of specific heats at equilibrium, c,; is heat capacity of it internal mode, c, is total 
heat capacity of the gas, 7; is the relaxation time of that internal mode, and the summation 
is performed over all internal degrees of freedoms (i.e., rotational, vibrational). However, 
the applicability of this expression is limited to the low frequency regime where wt << 1; w 
being the frequency of sound wave. [60]. 

Li et al. related bulk viscosity to bulk modulus and relaxation time as follows: 


Hb = K Tiot (10) 


where, K is bulk modulus of the fluid, defined as K = —V(OP/OV), and To is total average 
relaxation time of internal energy in all excited modes. 

Kustova et al. questioned a priory split of bulk viscosity in rotational and vibrational 
components as done by Tisza(3| and Cramer (21). They argued that this splitting lead to 
overprediction of bulk viscosity at low temperatures and bulk to shear viscosity ratio of 
~ 2000 for carbon dioxide at low temperatures is unjustified. They first derived a expression 
for bulk viscosity using modified Chapman-Enskog method under the assumption of 
local thermodynamic equilibrium 
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where Cint and Tini are total heat capacity of internal degrees of freedoms and corresponding 
relaxation time, respectively. Then following the work by Mason and Monchick[63], they 
assumed that collisions with simultaneous exchange of rotational and vibrational energies 


are rare. By invoking this assumption, it can be shown that 
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By using Eq and (12), bulk viscosity can be expressed as 
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On the basis of above formulation, Kustova et al. also deduced that bulk to shear viscosity 
ratio for CO, at 300-1000 K should be in the range 1 to 3. 


3.1.2 Experimental methods 


The experimental determination of bulk viscosity is not as straightforward as shear viscosity 
and usually based on indirect techniques, such as absorption and dispersion of sound wave, 
and Rayleigh-Brillouin scattering (65). 

Absorption of the sound wave is the additional decrease in intensity with distance, over 
and above the geometric reduction caused by the inverse square law. It has been found that 
the experimentally observed absorption is much higher than the predictions based on the 
theory that accounts only for classical absorption, i.e., absorption due to shear viscosity, ther- 
mal conductivity, and thermal radiation. Since this excess absorption cannot be attributed 
to dissipation phenomenon because of translational motion of molecules (i.e., shear viscosity, 
heat conduction), it is assumed that this excess absorption is because of bulk viscosity [66]. 
The absorption of sound is characterized by absorption coefficient (œ), and it is related to 
bulk viscosity, 4», as following [47]: 
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where, Peg is equilibrium pressure, M is molar mass, «& is thermal conductivity, and R is gas 
constant. 

However, the assumption that the excess absorption is due to bulk viscosity can only be 
examined when some direct measurements of bulk viscosity from an independent method 
are made and values are compared (66). Furthermore, this approach of measuring bulk 
viscosity is susceptible to considerable errors since it involves subtraction of classical absorp- 
tion coefficient from total absorption coefficient to get absorption due to bulk viscosity. For 
the calculation of classical absorption, the use of u and « taken from different sources also 
introduces error in estimates of bulk viscosity made using this method (65). 

Alternatively, bulk viscosity can also be measured by sound dispersion experiments. Dis- 


persion of sound causes speed of sound to be frequency dependent, and this dependency is 


given as follows (67|- 
(15) 


where a is the speed of sound at frequency w, ao and a. are speed of sound for very low 
and very high frequencies, respectively. The obtained relaxation time then can be related to 
bulk viscosity. 

Pan et al. suggested that bulk viscosity of dilute gases can also be measured using Co- 
herent Rayleigh — Brillouin Scattering (CRBS). In this technique, gas density perturbations 
are generated and measured using laser beams. The experimentally observed scattering 
profile is then compared with that obtained from the theoretical models to get transport 
coefficients, including bulk viscosity. However, in contrast to acoustic experiments, which 
measure bulk viscosity at megahertz frequencies, these experiments measure bulk viscosity in 
gigahertz frequencies. Because of this reason, a significant difference between the estimated 
values from the two above mentioned methods is usually observed [68H74]. 

Furthermore, Emanuel et al. |7| has deduced that for dense polyatomic gases, the density- 
based thickness of shock wave consists of many thousands of mean-free paths, and varies 
linearly with the ratio 4/u. Thus, ideally, the experimentally measured shock wave thickness 
can be used to calculate bulk viscosity; nevertheless, we could not find any experimental 


implementation of this technique in the literature. 


3.1.3 Computational methods 


The methods employed for the determination of transport properties using numerical sim- 
ulations are broadly classified in two categories- non-equilibrium simulation based meth- 
ods (e.g., non-equilibrium molecular dynamics (NEMD)), and equilibrium simulation based 
methods (e.g., equilibrium molecular dynamics (EMD)) (see Fig. [3). In the former ap- 
proach, the non-equilibrium responsible for the desired transport property is produced at 
the microscopic level, and then the transport property is related to other variables as in phys- 
ical experiments. The latter approach uses relations, such as Green-Kubo relation{75| [76], 
Einstein relation [78], or expressions derived from the Chapman-Enskog expansion[48} 


Eg 


Equilibrium based methods The Green-Kubo method uses the Green-Kubo relations 
for calculation of transport properties [r6]. For shear viscosity, u, and bulk viscosity, Hp, 
these relations are given as follows: 

y OO 


=T j, (Pal) PaO)dt (16) 


u 
where kg is the Boltzmann constant, T is temperature, P;;(t) denotes the instantaneous 


value of ijt? off-diagonal element of the pressure tensor at a time t, and the angle bracket 


indicates the ensemble average. Further, to reduce the statistical error in calculation of ju, 
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averaging is performed over three different values obtained from three different components 


of pressure tensor viz., Pij, Pjk, and Ppi. 


y oo 
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My = 
Here, P(t) is instantaneous value of the average of three diagonal terms of pressure tensor 
at a time t, i.e., P(t) = §[Py(t) + Pjj(t) + Pre(t)]. The fluctuations, 5P(t), is aberration of 
mean pressure from equilibrium pressure, i.e., dP(t) = P(t) — Peq; where Peq is equilibrium 
pressure of the system, and it is calculated by time average of P(t) over a long time. 

The Green-Kubo method is a robust way to measure transport coefficients, but some- 
times, they suffer from several issues. For example, for the correct estimation of viscosity 
coefficients, the auto-correlation function should decay to zero with time. In such a case, the 
integral of auto-correlation function would reach a constant value. Nevertheless, it does not 
necessarily happen in practice. The auto-correlation function might show either long-time 
tails or fluctuations [83]. Further, viscosities should be estimated from the region of the 
graph of the integral (of auto-correlation function) vs. time, when it reaches a constant 
value. First of all, it is difficult to identify such a region, and second, even if we can identify 
such a region, there will always be arbitrariness in the value of viscosity because of ambiguity 
in determining the cut-off time [83]. 

The Einstein relations also find their origin in linear response theory. These rela- 
tions relate shear and bulk viscosity to the slope of generalized mean-squared displacement 


functions as given below 
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Here, N is total number of particles, m is particle mass, r(t) and v, (t) are the position and 
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velocity vectors of n“ particle at time t, and 7;,;(t) and @,;(¢) are the Cartesian components 
of position and velocity vectors in direction 72 and j. 

The Einstein relations are theoretically equivalent to Green-Kubo relations. However, 
they can not be directly implemented in molecular dynamics (MD) simulations with periodic 
boundaries|86). The Einstein relations implicitly assumes the particles to follow a continuous 
trajectory. Whereas, in case of periodic boundaries, particles regularly exit from one face and 
enter from opposite face of the domain, therefore, the trajectory becomes discontinuous. This 


problem is overcome by using modified Einstein relations where the generalized displacement 
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functions are replaced with time integral of pressure tensor components. The modified 


Einstein relations are given as follow 
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The Chapman-Enskog theory based expressions for transport coefficients of dilute poly- 


atomic gases, derived by Taxman [81], requires the calculation of collision integrals. These 
collision integrals are usually evaluated through the Monte Carlo quadrature method[87], 
in which proper equilibrium distribution function is used to sample pre-collision state of 
molecules, and then post-collision states are calculated by solving molecular trajectory clas- 
sically [80]. The Chapman-Enskog method is computationally more efficient than the Green- 
Kubo method, however, the applicability of this method depends upon the availability of 
explicit expressions, whereas the Green-Kubo method, though less efficient, does not have 
any such constraint [50]. 

In another approach, transport coefficients can also be estimated by numerical simulation 
of the Rayleigh-Brillouin scattering experiment and its coherent version, coherent Rayleigh 
Brillouin scattering. In this method, a gas is simulated at equilibrium, and density fluctua- 
tions are sampled as time-series data. The discrete power spectrum is then obtained using 
the square of the Fourier transformation of this time series [89]. 

All the three above-mentioned methods, i.e., Chapman-Enskog expressions, Green-Kubo 
relations, and simulation of Rayleigh-Brillouin scattering can be employed in either MD and 


classical trajectory direct simulation Monte Carlo (CT-DSMC) simulation framework. More 


details on these three methods are available in Ref. [74] (90). 


Non-equilibrium based methods Non-equilibrium based methods measures transport 
coefficients by directly measuring the gradient of the corresponding parameter. Such methods 
are well developed for shear viscosity, heat conductivity, and mass diffusivity. However, for 
bulk viscosity, historically it has been seen that implementation of such a method is a 
challenging task. To the best of our knowledge, only one attempt has been made so far to 
use non-equilibrium based methods for the estimation of bulk viscosity. In this work, Hoover 
et al. cyclically compressed and expanded the fluid in molecular dynamics framework 


in the following manner to produce measurable effects: 


L/Lo = 1+ €sin(wt) (22) 
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where, Lo is mean length, L is instantaneous length of the cubic simulation domain, € is 


strain amplitude, and w is the frequency describing the linearized strain rate. 
é = fw cos(wt) (23) 


As the linearized strain rate (€) approaches zero, authors expected the average pressure of 
the system to deviate from equilibrium pressure by —3&wyip cos(wt). Also, if the deformation 
given by Eq. takes place through the mechanism of external work, then the lost work 
due to irreversible heating will give rise to energy increase per cycle, AE per cycle, (the cycle 


time being 27/w), which is equal to 
AFEper cycle = (2r /w)9E’ w uV /2 (24) 


Based on the measurement of average pressure and rise in the energy of the system, 
authors estimated bulk viscosity of a soft sphere fluid modeled with potential ¢(r) = e(o /r)!?. 

All of the methods discussed above have one or more limitations. The sound absorption 
and dispersion method give frequency-dependent bulk viscosity values. However, since bulk 
viscosity is a transport coefficient, it should depend only on the state of fluid, not on the pro- 
cess it is undergoing. The optical method, i.e., Rayleigh-Brillouin scattering, can not account 
for vibrational effects, since the frequency used is of GHz range, and 1/frequency becomes 
smaller than vibrational relaxation time. In the second vertical, i.e., numerical methods, 
the Green-Kubo method sometimes faces difficulties in convergence due to long-time tails. 
Chapman-Enskog method is very robust and fastest among all numerical methods. However, 
Chapman-Enskog relations are difficult to obtain except for simple cases like monatomic and 
diatomic gases. In the second sub-vertical, there is only one nonequilibrium method, i.e., 
cyclic compression and expansion by Hoover et al.[19}. There is only one application avail- 
able in the literature[9 1]. Since this method uses energy absorbed over many compression 
and expansion cycles to estimate bulk viscosity, it may also give frequency-dependent bulk 
viscosity values. Moreover, this method does not give much insights associated with physics 
of one single cycle. Therefore, to address these issue, Sharma et al. proposed a first 
principle based continuous compression/expansion method which can determine bulk vis- 
cosity directly from the difference between mechanical (Pmecn) and thermodynamic pressure 
(Pihermo): 

Their method uses numerical measurement of Pmecen and Pihermo in an NEMD simu- 
lation of an expanding fluid, and then relates the bulk viscosity to them by the relation 
ua = (Pihermo — Pmecnh)/V - U, where V - Ñ is the controlled rate of expansion of the fluid 
per unit volume. The key success of this method was that it, being inherently based on 


expansion /compression of the gas, allowed them to investigate the effects of nonequilibrium 
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on bulk viscosity, e.g., variation of bulk viscosity with magnitude (|V - w|) and direction of 
volumetric change (i.e., expansion vs. compression). Hence, this method enabled them to 


gain detailed insights of associated flow physics. 


3.2 Determination of bulk biscosity of QGP and hadronic matter 


For the estimation of transport coefficient of QGP and hadronic matter, two standard ap- 
proaches are mainly used — the Boltzmann equation based relaxation time approximation 
(RTA) approach and the linear response theory based Kubo/Green-Kubo formulation. A 
brief review of these methods can be found in Ref. [52]. A vast amount of research 
has been done on this topic. Here we summarize only some of the key contributions in this 
field. Gavin used the well known non-relativistic form of the Boltzmann equation to 
calculate the transport coefficients for both the QGP and hadronic matter using the relax- 
ation time approximation (RTA) method. Prakash et al. studied the equilibration of 
hot hadronic matter in the framework of relativistic kinetic theory. They calculated trans- 
port coefficients considering only elastic collisions in the dilute gas limit using extended 
Chapman-Enskog formalism. For a general review of the relativistic kinetic theory, reader 
may refer to the classical text by Groot et al. [56]. Chakraborty et al. extended 
the classical works of Prakash et al. and Gavin[53], and presented a theoretical 
framework for the calculation of shear and bulk viscosity of hot hadronic matter. Their 
work accounted for not only inelastic collisions but also formation and decay of resonances, 
temperature-dependent mean fields, and temperature-dependent masses. Demir et al. 
carried out Ultra-relativistic Quantum Molecular Dynamics (UrQMD) simulations of 
hadronic media, and calculated bulk viscosity using both Green-Kubo method and relaxation 


time approximation. 


4 Applications of bulk viscosity 


To evaluate the validity of Stokes’ hypothesis, let us consider a supersonic flow of air and 
analyze the terms of the equation Pyech = Pihermo — Hb V - U. For bulk viscosity effects to 
be significant, uV - t should be comparable to Pinermo. Assuming pressure is ~ 10° Pa and 
L» is ~ 107% Pa s, even if V - Ẹ is as high as © 104 s™t, the difference between mechanical 
and thermodynamic pressure would be just 0.1 Pa. This difference Pinermo — Pmech = 0.1 
Pa is six orders of magnitude smaller than the Pipermo. Therefore, in most of the commonly 
encountered flow problems, it is safe to assume bulk viscosity to be zero. However, bulk 
viscosity effects may become important when either the V - w is very high (e.g., inside 


a shock wave), or when fluid is compressed and expanded in repeated cycles such that the 
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cumulative effect of the small contributions from each cycle is no more negligible (e.g., sound 
wave) [22], or when the atmosphere consists of the majority of those gases, such as CO», 
which exhibit a large bulk viscosity 5l, or when results of interest might get affected by 
even small disturbances, e.g., the study of Rayleigh-Taylor instability [15]. In such cases, it 
becomes necessary to account for the bulk viscosity terms in the Navier-Stokes equation. 

Several researchers have investigated the effects of the incorporation of bulk viscosity in 
analytical or CFD studies of various flow scenarios. Emanuel et al. reviewed 
bulk viscosity and suggested that the effects of bulk viscosity should be accounted for in the 
study of high-speed entry into planetary atmospheres. They observed that the inclusion of 
bulk viscosity could significantly increase heat transfer in the hypersonic boundary layer[5). 
Chikitkin studied the effects of bulk viscosity in flow past a spacecraft. They reported 
that the consideration of bulk viscosity improved the agreement of velocity profile and shock 
wave thickness with experiments. Shevlev studied the effects of bulk viscosity on CO%2 
hypersonic flow around blunt bodies. The conclusions of their study were in line with that 
of Emanuel. They suggested that incorporation of bulk viscosity may improve predictions 
of surface heat transfer and other flow properties in shock layer. 

Elizarova et al. and Claycomb et al. carried out CFD simulations of normal 
shock. They found that including bulk viscosity improved the agreement with experimental 
observations for shock wave thickness. A recent study by Kosuge and Aoki on shock- 
wave structure for polyatomic gases also confirms the same. Bahmani et al. studied 
the effects of large bulk to shear viscosity ratio on shock boundary layer interaction. They 
found that a sufficiently high bulk to shear viscosity ratio can suppress the shock-induced 
flow separation. Singh and Myong studied the effects of bulk viscosity in shock-vortex 
interaction in monatomic and diatomic gases. They reported a substantially strengthened 
enstrophy evolution in the case of diatomic gas flow. Singh et al. investigated the 
impact of bulk viscosity on the flow morphology of a shock-accelerated cylindrical light bub- 
ble in diatomic and polyatomic gases. They found that the diatomic and polyatomic gases 
have significantly different flow morphology than monatomic gases. They produce larger 
rolled-up vortex chains, various inward jet formations, and large mixing zones with strong, 
large-scale expansion. Touber studied the effects of bulk viscosity in the dissipation of 
energy in turbulent flows. He found that large bulk-to-shear viscosity ratios may enhance 
transfers to small-scale solenoidal kinetic energy, and therefore, faster dissipation rates. Ri- 
abov questioned the ability of bulk viscosity to model spherically-expanding nitrogen 
flows in temperature range 10 to 1000 K by comparing results to Navier-Stokes equations 
to relaxation equation. He reported that the bulk viscosity approach predicts much thinner 
spherical shock wave areas than those predicted by relaxation equations. Moreover, the dis- 


tributions of rotational temperature along the radial direction predicted by the bulk viscosity 
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approach had neither any physical meaning nor matches with any known experimental data 
for expanding nitrogen flows. 

Fru et al. performed direct numerical simulations (DNS) study of high turbulence 
combustion of premixed methane gas. They found that the incorporation of bulk viscosity 
does not impact flame structures in both laminar and turbulent flow regimes. Later, the same 
group extended their study to other fuels, viz., hydrogen and synthetic gas. In this study 
[11], they found that though flame structures of methane remained unchanged before and 
after incorporation of bulk viscosity, the same for hydrogen and syngas showed noticeable 
modifications. Sengupta et al. studied the role of bulk viscosity on Rayleigh Taylor 
instability. They found that the growth of the mixing layer depends upon bulk viscosity. 
Pan et al. has shown that bulk viscosity effects cannot be neglected for turbulent flows 
of fluids with high bulk to shear viscosity ratio. They found that bulk viscosity increases 
the decay rate of turbulent kinetic energy. Boukharfane et al. studied the mechanism 
through which bulk viscosity affects the turbulent flow. They found that the local and 
instantaneous structure of the mixing layer may vary significantly if bulk viscosity effects 
are taken into account. They identified that the mean statistical quantities, e.g., the vorticity 
thickness growth rate, do not get affected by bulk viscosity. On this basis of their study, 
they concluded that results of refined large-eddy simulations (LES) might show dependence 
on the presence/absence of bulk viscosity, but Reynolds-averaged Navier-Stokes (RANS) 
simulations might not as they are based on statistical averages. 

Connor studied the effects of bulk viscosity in the compressible turbulent one-, two-, 
and three-dimensional Couette flows through DNS simulations. The objective of the study 
was to test whether invoking the Stokes’ hypothesis introduces significant errors in the anal- 
ysis of compressible flow of solar thermal power plants, and carbon capture and storage 
(CCS) compressors. They found that most of the energy is contained in the solenoidal veloc- 
ity for both CCS and concentrated solar power plants. Therefore, assuming bulk viscosity 
to be zero does not produce any significant errors, despite the compressors operating at 
supersonic conditions. However, bulk viscosity effects may become significantly close to the 
thermodynamic critical point. 

Billet et al. showed that the inclusion of bulk viscosity in CFD simulations of 
supersonic combustion modifies the vorticity of the flow. Lin et al. have shown 
that acoustic wave attenuation in CFD simulations can be accounted for by incorporating 
bulk viscosity. Nazari studied the influence of liquid bulk viscosity on the dynamics of a 
single cavitation bubble. They reported that bulk viscosity significantly affects the collapse 
phase of the bubble at high ultrasonic amplitudes and high viscosities. High bulk viscosity 
values also altered the maximum pressure value inside the bubble. 


The classical Navier-Stokes equation loses its accuracy as the deviation from equilibrium 
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increases, or the extent of nonequilibrium is high. This usually happens in the study of 
high-temperature gas dynamics and rarefied gas dynamics. There are primarily three mod- 
els used in the continuum framework to model nonequilibrium, viz., state-to-state, multi- 
temperature, and one-temperature model. The state-to-state model is computationally most 
costly and can be used in systems that are far from equilibrium and have coupled fluid, ther- 
mal and chemical kinetics. In comparison, the one-temperature model is computationally 
most simple and suitable for near-equilibrium systems. In one temperature model, both ro- 
tational and vibrational relaxations are accounted for using bulk viscosity coefficient. On the 
other hand, in the state-to-state approach, the vibrational chemical kinetics is described by 
master equations for populations of vibrational states, and the fast rotational relaxation is 
accounted for through bulk viscosity parameters. Similarly, in the multi-temperature model, 
vibrational kinetics is governed by relaxation equations for various vibrational modes, and 
the rotational relaxation is modeled through bulk viscosity. For a more detailed discussion 
on these models, the reader is referred to the Ref. [38]. 

In addition to these classical fluid-dynamics scenarios, bulk viscosity may also play an 
important role in several cosmological phenomena, e.g., damping of vibrations created during 
the formation of a new neutron star, and growth of gravitational—wave instability in rapidly 
rotating neutron star (40). Origin of bulk viscosity in these circumstances is primarily 
due to the chemical nonequilibrium caused by nuclear reactions. Bulk viscosity along with 
other transport properties is also of central importance to the space-time description of the 
heavy-ion collision experiments being conducted at the Brookhaven National Laboratory’s 
Relativistic Heavy Ion Collider (RHIC) and CERN’s Large Hadron Collider (LHC). One of 
the primary objectives of these experiments was the formation and investigation of Quark- 


Gluon-Plasma (QGP), the state of the matter of the Early Universe (first 30 ys after the 
Big Bang). 


5 Concluding remarks and future scope 


Although of great importance in several fluid mechanical phenomena, bulk viscosity is still 
one of the gray areas of fluid mechanics. Despite the fact that numerous theoretical, exper- 
imental, and computational studies have been carried out in the past, there is still lack of 
well-established values for bulk viscosity of common gases such as nitrogen, oxygen, and car- 
bon dioxide. For most gas mixtures, data is not available at all, or if available, experimental 
or numerical, is spread in a wide range. Therefore, studies that determine the bulk viscosity 


of common fluids are needed to achieve a common consensus on accepted values. 
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Disclaimer 


The contents of this article is primarily based on Introduction chapter of the thesis: Bhanu- 
day Sharma. “On the Nature of Bulk Viscosity of Dilute Gases”. PhD thesis. Indian 
Institute of Technology Kanpur, 2022 


Appendix A: Selected references for further reading 


Reader is redirected to following resources for further understanding of concepts and his- 
toric developments in the subject of bulk viscosity. Enumerated items are arranged in the 


suggested order of reading. 
e For first introduction to the subject: 
1. Mohamed Gad-el Hak. “Questions in fluid mechanics: Stokes’ hypothesis for a 
Newtonian, isotropic fluid”. In: Journal of Fluids Engineering 117.3 (1995), p. 5 
2. Guido Buresti. “A note on Stokes’ hypothesis”. In: Acta Mechanica 226.10 
(2015), pp. 3555-3559 
e For a detailed understanding: 
1. Bhanuday Sharma. “On the Nature of Bulk Viscosity of Dilute Gases”. PhD 
thesis. Indian Institute of Technology Kanpur, 2022 


2. Lev Davidovich Landau and Evgenii Mikhailovich Lifshitz. Course of Theoretical 
Physics - Fluid Mechanics. Oxford, England: Elsevier, 2013 (Chapter 2, Sec 15 
and Chapter-8, Sec 81) 


3. Karl F Herzfeld and Theodore A Litovitz. Absorption and dispersion of ultrasonic 
waves. Vol. 7. Academic Press, 2013 


e For historic developments: 


1. Leonard Hall. “The origin of ultrasonic absorption in water”. In: Physical Review 
73.7 (1948), p. 775 


2. SM Karim and L Rosenhead. “The second coefficient of viscosity of liquids and 
gases”. In: Reviews of Modern Physics 24.2 (1952), p. 108 


3. L Rosenhead. “A discussion on the first and second viscosities of fluids”. In: 
Proc. R. Soc. London 226 (1954), pp. 1-69 
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. RO Davies. “Kinetic and thermodynamic aspects of the second coefficient of vis- 


cosity”. In: Proceedings of the Royal Society of London. Series A. Mathematical 
and Physical Sciences 226.1164 (1954), pp. 24-34 


. SV Borovinskii and OV Mazurin. “On the concept of volume viscosity”. In: 


Progress and Trends in Rheology II. Springer, 1988, pp. 79-81 


. C Truesdell. “The present status of the controversy regarding the bulk viscosity 


of fluids”. In: Proc. R. Soc. Lond. A 226.1164 (1954), pp. 59-65 


. George Emanuel. “Bulk viscosity of a dilute polyatomic gas”. In: Physics of 


Fluids A: Fluid Dynamics 2.12 (1990), pp. 2252-2254 


. Willard E Meador, Gilda A Miner, and Lawrence W Townsend. “Bulk viscosity 


as a relaxation parameter: fact or fiction?” In: Physics of fluids 8.1 (1996), 
pp. 258-261 


. Hui Dong, Nan Su, and Qun Wang. “Bulk viscosity in nuclear and quark matter”. 


In: Journal of Physics G: Nuclear and Particle Physics 34.8 (2007), 5643 


e For complete treatment of equation of motion of viscous fluids without Stokes’ hypoth- 


esis: 


1. 


George Emanuel. Analytical fluid dynamics. CRC press, 2000 


e For fundamentals of equilibrium statistical mechanics: 


1. 


Biman Bagchi. Statistical Mechanics for Chemistry and Materials Science. CRC 
Press, 2018 


. Michael P Allen and Dominic J Tildesley. Computer simulation of liquids. Oxford 


university press, 2017 
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